#!/usr/bin/env python
import pyferret

infl = "./Data/Fig10b.nc"

pyferret.start(quiet=True)

pyferret.run("set window/aspectt=.9")

pyferret.run("use \"%s\""  % (infl))
pyferret.run("set region/l=32:44")


pyferret.run('fill/nolab/nokey/level=(0,12,1)(+inf) precip')
pyferret.run('ppl shakey 1, 0, .12, 0, 3, 9,`($ppl$xlen)-7.0`,`($ppl$xlen)-7.0+8.25`,`($ppl$ylen)-5.7`,`($ppl$ylen)-5.7+0.15`')
pyferret.run('ppl fill')

pyferret.run('contour/o/nolab/level=(-0.5,0.5,0.1) slope*10')
pyferret.run('contour/o/nolab/color=red/size=0.00001/pen=8/level=(0.2) sig')

title="(b) HiFLOR (AllForc, Lon=125-142E, 1979-2015)"
pyferret.run('label/nouser `($ppl$xlen)/2`, `($ppl$ylen) + 0.2`,0,0,0.20, %s' % (title))

pyferret.run('label/nouser `($ppl$xlen)/2.0+0.15`,`($ppl$ylen)-6.85` 1, 0, 0.14, Date')
pyferret.run('label/nouser `($ppl$xlen)/2.0+2.0`,`($ppl$ylen)-7.5` 1, 0, 0.12, Climatological mean precipitation (mm/day)')

pyferret.run('label/nouser `($ppl$xlen)-8.55`,`($ppl$ylen)/2` 1, 90, 0.14, Latitude')

pyferret.run("frame /file=\"Fig10b.png\"")
pyferret.run("quit")
